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Abstract. Nuclear reaction rates are quantities of fundamental importance in astrophysics. Substantial efforts 
have been devoted in the last decades to measure or calculate them. The present paper presents for the first 
time a detailed description of the Brussels nuclear reaction rate library BRUSLIB and of the nuclear network 
^ • generator NETGEN in a journal and in a format so as to make these nuclear data packages easily accessible 

to astrophysicists for a large variety of applications. BRUSLIB is made of two parts. The first one contains 
the 1999 NACRE compilation based on experimental data for 86 reactions with (mainly) stable targets up to 
Si. BRUSLIB provides an electronic link to the published, as well as to a large body of unpublished, NACRE 
data containing adopted rates, as well as lower and upper limits to these rates. The second part of BRUSLIB 
concerns nuclear reaction rate predictions that have to complement the experimentally-based rates for use in a 
large variety of astrophysics modellings. An electronic access is provided to tables of thousand of rates calculated 
within a statistical Hauser-Feshbach approximation, which limits the reliability of the rates to reactions producing 
compound nuclei with a high enough level density. These calculations make use of global and coherent microscopic 
Oh' nuclear models for the quantities entering the rate calculations. The use of such models is utterly important, 

and makes the BRUSLIB rate library unique. A description of the Nuclear Network Generator NETGEN that 
5_l ■ complements the BRUSLIB package is also presented. NETGEN is a tool to generate nuclear reaction rates for 

temperature grids specified by the user. The information it provides can be used for a large variety of applications, 
including Big Bang nucleosynthesis, the energy generation and nucleosynthesis associated with the non-explosive 
and explosive hydrogen to silicon burning stages, or the synthesis of the heavy nuclides through the s-, a- and r-, 
rp- or p-processes. 
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1. Introduction 

Since around the nineteen fifties, astrophysics has ad- 
vanced at a remarkable pace, and has achieved an impres- 
sive record of success. One of the factors contributing to 
these rapid developments is undoubtedly a series of spec- 
tacular breakthroughs in nuclear astrophysics, which em- 
bodies the special interplay between nuclear physics and 
astrophysics. 

The close relationship between these two major scien- 
tific disciplines comes about because of the clear demon- 
stration that the Universe is pervaded with nuclear physics 
imprints at all scales. Starting with the Big Bang nucle- 
osynthesis episode, the structure, evolution and composi- 
tion of a large variety of cosmic objects, including the Solar 
System and its various constituants (down to meteoritic 
grains), bear strong imprints of the properties of atomic 
nuclei, as well as of their interactions. Therefore, careful 
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and dedicated experimental and theoretical studies of a 
large variety of nuclear processes are indispensable tools 
for the modeling of the ultra-macroscopic systems astro- 
physics has to deal with. 

Over the years, an impressive body of nuclear data of 
astrophysics interest have been obtained through labora- 
tory efforts, complemented with theoretical developments. 
The latter are indispensable as nuclear experiments con- 
ducted today or expected to be performed in any foresee- 
able future can and will cover only a minute fraction of 
the astrophysics needs. This is because an extremely large 
variety of highly unstable 'exotic' nuclei that, for long, 
will not be produced and studied in the laboratory are 
expected to be involved in the modelling of many astro- 
physics processes and events. Many of the basic properties 
of these nuclei are to be known for this purpose, and so 
are their interactions, in particular those with nucleons or 
a-particles. Even when laboratory-studied nuclei are con- 
sidered, theory has very often to be called for assistance. 
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In many respects, laboratory conditions are indeed very 
different from stellar ones, which are highly versatile and 
are often charaterized by high temperatures and/or densi- 
ties that are out of reach of laboratory simulations. In ad- 
dition, nuclear reactions between charged particles inside 
non-exploding stars take place in an energy regime that is 
in all but a few exceptional cases out of reach of direct ex- 
perimental scrutiny. Indirect methods bring some comple- 
ment of information, but clearly do not cover all the needs. 
In explosive situations, the energies of astrophysical inter- 
est are higher, and the cross sections are correspondingly 
larger. However, in such events, many reactions involve 
unstable targets, so that the fraction of those reactions 
of potential interest for which experimental reaction data 
are lacking is even larger. 

The experimental and theoretical nuclear physics 
achievement mentioned above in fact let nuclear astro- 
physics face a new and difficult challenge. The rapidly 
growing wealth of nuclear data becomes, ironically, less 
and less easily accessible to the astrophysics community. 
Mastering this volume of information and making it avail- 
able in an accurate and usable form for incorporation into 
astrophysics models become urgent goals of prime neces- 
sity. The establishment of the required level in the priv- 
ileged communication between nuclear physicists and as- 
trophysicists necessarily requires the build-up of well doc- 
umented and evaluated sets of experimental data or theo- 
retical predictions of astrophysical relevance. Some years 
ago, the Institut d'Astronomie et d'Astrophysique of the 
Universite Libre de Bruxelles has decided to tackle this 
challenge. This work has grown into the BRUSsels nu- 
clear LIBrary for astrophysics applications, referred to as 
BRUSLIB, and into a nuclear NETwork GENerator called 
NETGEN. 

One goal of this paper is to present for the first time 
a detailed description of BRUSLIB and of NETGEN in a 
journal and in a format such that these nuclear packages 
become easily accessible to astrophysicists for a large vari- 
ety of applications. In particular, the BRUSLIB reaction 
rates presented in a tabular form that is well suited to as- 
trophysics needs are available electronically at the address 
http://www-astro.ulb.ac.be 

BRUSLIB is composed of two main parts. The first 
one concerns the Nuclear Astrophysics Compilation of 
REaction rates referred to as NACRE. This compilation 
provides the rates of 86 thermonuclear reactions of astro- 
physical relevance based on an in-depth analysis of exper- 
imental data. Its description is presented in Sect. The 
second part describes theoretical evaluations of a collec- 
tion of about 100000 rates for thermonuclear reactions in- 
duced by nucleons or a-particles, as well as photo-induced 
reactions not included in NACRE, including nuclei with 
8 < Z < 110 lying between the proton and the neutron 
drip lines (Sect.|3J). The rate calculations are based on the 
statistical Hauser-Feshbach (hereafter HF) model. They 
require the knowledge of a substantial amount of data con- 
cerning basic properties of the nuclei and of their interac- 
tions. The predictions of these properties rely on the use of 



global and coherent microscopic nuclear models. These ut- 
terly important model characteristics make the BRUSLIB 
rate library unique. A confrontation between selected ex- 
perimental data and the predictions used in BRUSLIB is 
provided in Sect.rj] 

A second goal is to describe NETGEN (Sect. EJ), a 
package for constructing nuclear reaction networks on 
grounds of the nuclear physics input from BRUSLIB and, 
when necessary, from other sources. 

It has to be noted that no information is provided 
by the present release of BRUSLIB on the rates of non- 
statistical ('direct') or of non-thermonuclear ('spallation') 
reactions that can develop in low-temperature and low- 
density astrophysical media, like the interstellar or cir- 
cumstellar medium. 

2. The experimentally-based NACRE evaluation 
and compilation of nuclear reaction rates 

The necessity of placing relevant nuclear data of high qual- 
ity at the disposal of the astrophysics community has been 
the driving motivation for the construction of a compila- 
tion aimed at superseding the work of Fowler and col- 
laborators (see Caughlan & Fowler 1988, hereafter CF88, 
and references therein to former compilations). The goal 
was not just to update the CF88 rates with newly avail- 
able experimental data, but also to modify quite deeply 
different aspects of the format of the CF88 compilation. 
At this point, slightly more than half of the CF88 rates 
have been re-compiled on grounds of a careful evalua- 
tion of experimental data available by 15 June 1998 , 
and the results of the work make the so-called NACRE 
(Nuclear Astrophysics Compilation of REactions) compi- 
lation (Angulo ct al. 1999). The reactions analyzed up to 
now are listed in Table 1. They comprise an ensemble of 
86 charged particle induced reactions on stable targets up 
to Si involved in Big Bang nucleosynthesis and in the non- 
explosive H- and He-burning modes, complemented with 
a restricted number of reactions of special astrophysical 
significance on the unstable 7 Be, 22 Na and 26 A1 nuclides. 
An updated and enlarged version of NACRE is currently 
in preparation. 

NACRE is described in detail by Angulo et al. (1999). 
It contains in particular 

(1) the formalism that has been adopted in order to de- 
rive the Maxwellian-averaged astrophysical rates of the 
(exothermic and endothermic) charged particle induced 
reactions and of their reverse; 

(2) a general description of the treatment of the data, 
which has sometimes to be adapted to specific cases, as 
described in the comments accompanying each of the com- 
piled reactions. If necessary, tables of narrow resonances 
with their characteristics are provided. A careful analy- 
sis of the experimental uncertainties in the quantities in- 
volved in the evaluation of the reaction rates is carried out 
for each reaction. From this, adopted rates are provided, 
as well as low and high limits. This large-scale estimate of 
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Table 1. List of the NACRE compiled reactions. References to compilations, new data or evaluations concerning some 
of these reactions that have appeared after the completion of NACRE are provided. A + sign indicates that only the 
most recent reference to a given reaction is provided. Reaction rates are given in the corresponding reference (a) in 
tabular form, (b) as analytic formulae, or (c) as a figure. 
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the uncertainties is considered to be an essential feature 
of the NACRE compilation; 

(3) the procedure adopted for extrapolating experimen- 
tal data when required in order to evaluate the reaction 
rates. This extrapolation to the very low energies neces- 
sary to evaluate the rates down to the lowest considered 
temperatures raises many difficult problems. A particu- 
lar one concerns the necessity to correct the laboratory 
cross sections for electron screening (e.g. Rolfs & Rodney 



1988) before their use in the rate calculations. This ef- 
fect becomes significant for relative energies E of the re- 
action partners such that E/U e < 100, where U e is the 
so-called screening potential. An approximate procedure 
is devised in order to eliminate data that can be 'polluted' 
by laboratory screening. On the other hand, in stellar con- 
ditions, the nuclei are surrounded by a dense electron gas 
that reduces the Coulomb repulsion and makes the pene- 
tration of the Coulomb barrier easier. The cross sections 
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are therefore enhanced in comparison with those of reac- 
tions between bare nuclei. This stellar screening effect can 
be evaluated by applying, for example, the Debye-Hiickel 
theory (e.g. Cox & Giuli 1968). The NACRE rates exclude 
these stellar screening factors. 

The evaluation of the rates for temperatures as high as 
Tg = 10 also requires the application of special extrapola- 
tion techniques when reliable cross section data are lacking 
at sufficiently high energies. In those cases, a duly docu- 
mented HF approach is used (see also Sect. 13, 1|) - which 
smoothly connects to the experimentally-based rate esti- 
mates; 

(3) the procedure used to evaluate the contribution of ex- 
cited states of the target nuclei to the effective stellar re- 
action rates. In a stellar plasma, the excited levels of a tar- 
get nucleus are indeed thermally populated, and thus con- 
tribute to the reaction mechanism. As a result, the stellar 
rates may differ from those obtained when the target nu- 
clei are in their ground state. This difference is expressed 
in terms of a correction factor ru that has to multiply the 
target ground state rate in order to obtain the stellar rate. 
In general, this correction cannot be derived experimen- 
tally. The very approximate treatment of this correction 
by Caughlan & Fowler (1988) is replaced in NACRE by a 
more quantitative procedure based on the use of the HF 
model (Sect. I3.1|) . and on the classical assumption of a 
Maxwellian population of the nuclear excited states. Note 
that this assumption may be invalid if the target nucleus 
has an isomeric state which is not in thermal equilibrium 
with the ground state in certain astrophysical situations. 
Such an example of astrophysical interest concerns 26 A1, 
to which a special treatment is applied. 

The calculation of r tt requires in particular the knowl- 
edge of the temperature-dependent partition functions of 
the target nuclei normalized to their ground state 

^) = ^£(2Jf + l)exp(-|), (1) 

where the summation extends over the states y, of target 
i with excitation energy and spin-parity J? (/i = for 
the ground state); 

(4) The values of the adopted, low and high Maxwellian- 
averaged unscreened rates on the ground state target are 
provided in tabular form for a selection of temperatures 
in the 0.001 < Tg < 10 range (Tg is the temperature in 
billion K). As they may be desired by some users, analyti- 
cal formulae approximating the tabulated adopted ground 
state rates or those including the contribution of thermal- 
ized excited states are provided as well. It has to be noted 
that some of these expressions have a form deviating from 
those classically advertised in the nuclear astrophysics lit- 
erature (e.g. Caughlan & Fowler 1988; Rolfs & Rodney 
1988; see e.g. Eq. (23) of Angulo et al. 1999). They are 
considered to provide more secure approximations to the 
numerically calculated Maxwellian-averaged rates. In ad- 
dition, NACRE provides analytical formulae for the mul- 
tiplicative factor that has to be applied to the rate of a 



given reaction in order to evaluate the rate of the inverse 
process, as well as for the nuclear partition functions for 
the targets involved in the compiled reactions (see Eq. ^) . 

NACRE is accessible electronically through the 
BRUSLIB library website http://www-astro.ulb.ac.be. 
Much more material is available there than presented 
above or published by Angulo et al. (1999). In particular, 
reaction rates are available for far more extended temper- 
ature grids. In addition, graphical material (astrophysical 
S-factors) can be retrieved. 

Finally, it is noted that other astrophysics-oriented 
experimentally-based reaction rate compilations have ap- 
peared around the time of the NACRE publication 
(Adelberger et al. 1998), or later (Iliadis et al. 2001; 
Descouvemont et al. 2004). The reactions considered in 
the latter two compilations are identified in Table 1, which 
also provides references to additional data or evaluations 
of relevance. This additional information will of course 
be duly analyzed and integrated into the new version of 
NACRE currently in preparation. 

3. The BRUSLIB theoretical reaction rates 

Much effort has been devoted in the last decades to mea- 
surements of reaction cross sections for astrophysical pur- 
poses. As already stated in Sect. 2] difficulties related to 
the conditions prevailing in the astrophysical plasmas re- 
main, however. In particular, charged-particle induced re- 
actions at stellar energies (far below the Coulomb barrier) 
have extremely small cross sections that are highly diffi- 
cult to measure. Specific difficulties are also raised by the 
measurements of photoreactions that have to be included 
in many astrophysics models (e.g. Arnould & Goriely 2003 
for a review). In addition, a huge number (thousands) of 
reactions of relevance involve more or less exotic nuclides. 
Clearly, many of these reaction rate experiments will re- 
main unfeasible for a long time to come. Theory has thus 
to supply the necessary data, which represents a major 
challenge of its own. 

BRUSLIB provides astrophysicists with a very ex- 
tended set of thermonuclear rates in the 0.01 < Tg < 10 
temperature range for all the reactions induced by neu- 
tron, proton and cv-particle captures by all nuclei with 
Z, N > 8 and Z < 110 located between the neutron and 
proton drip lines (i.e some 8000 nuclei). The rates of the 
(7,n), (7,p) and (7,0;) photodisintegrations of all these nu- 
clides are also tabulated for the same temperature grid. 
The calculations rely on the code MOST (Goriely 1998; 
Arnould & Goriely 2003) based on the HF model (Hauser 
& Feshbach 1952). This approach is valid if the density 
of nuclear levels of the compound systems formed as a re- 
sult of these captures is high enough. This is indeed the 
case if the targets are heavy enough and if they are lo- 
cated far enough from the proton or neutron drip lines to 
ensure that the excitation energies of the compound sys- 
tems are high enough. Roughly speaking, the reliability 
of the BRUSLIB rates may thus be limited to nuclides 
with mass numbers A > 40 close enough to the valley of 
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stability, this mass limit being of course shifted to higher 
and higher values as one moves further and further away 
from the valley. If these constraints are not met, a non- 
statistical treatment is more suited (Goriely 1997). It has 
to be noticed that purely theoretical HF rates are provided 
for some reactions already included in NACRE (Table 1). 
In these cases, it is advised to adopt the NACRE rates, or 
those from the Iliadis et al. (2001) or Descouvemont et al. 
(2004) compilations, as these rely on experimental data. 

A few words are in order at this point to clarify 
the basic philosophy governing the selection of the mod- 
els necessary for the calculation of the various ingredi- 
ents entering the HF BRUSLIB rates. Their evaluation is 
based on global and coherent microscopic (or at least semi- 
microscopic) models. These features have an importance 
in the astrophysics modellings that cannot be underesti- 
mated, and make the BRUSLIB rate library unique. They 
are indeed not shared by HF calculations based on other 
codes, like Talys (Koning et al. 2002), Empire (Herman et 
al. 2002), or Non-Smoker (Rauscher & Thielemann 2001). 

The global character of the underlying models is made 
highly desirable by the fact that, for many specific appli- 
cations including nuclear astrophysics, a very large body 
of reaction rates for which no experimental data exist have 
to be provided. The microscopic nature of the underlying 
models is essential as well. For nuclear astrophysics, as 
well as in other fields, a large amount of data need to be 
extrapolated far away from experimentally known regions. 
In these situations, two major features of the nuclear the- 
ories have to be considered. The first one is the accuracy 
of a model. In most nuclear applications, this criterion 
has been the main, if not the unique, one for selecting a 
model. The second one is the reliability of the predictions. 
A physically sound model that is based on first principles 
and is as close as possible to a microscopic description of 
the nuclear systems is expected to provide the best possi- 
ble reliability of extrapolations far away from experimen- 
tally known regions. Of course, the accuracy of such micro- 
scopic models in reproducing experimental data may be 
poorer than the one obtained from more phenomenologi- 
cal models in which enough free parameters can guaran- 
tee a satisfactory reproduction of the data at the expense 
of the quality of the input physics, and consequently of 
the reliability. The coherence (or 'universality') of these 
microscopic models (through e.g. the use of the same ba- 
sic nuclear inputs, like the effective nuclear forces) is also 
required as different ingredients have to be predicted in 
order to evaluate each reaction rate. Failure to meet this 
requirement could lead to bad ways in the rate evalua- 
tions. 

The BRUSLIB MOST rates are based on microscopic 
(or at least semi-microscopic) models for the necessary 
nuclear structure ingredients or interaction properties 
(Sect. I3.1|) that have been constructed to be global and 
coherent to the largest possible extent. In addition, these 
models have the major advantage of reaching a satisfac- 
tory compromise between accuracy and reliability for the 
relevant nuclear inputs, and consequently for the rates 



themselves. The level of accuracy of these models is in fact 
fully comparable to the one obtained from available more 
or less highly parametrized phcnomcnological approaches, 
while their reliability is by far better. 

Finally, the BRUSLIB reaction rates duly take into ac- 
count the necessary astrophysical specificities. The tem- 
perature dependence of the rates is predicted from the 
consideration of the Maxwell-Boltzmann distribution of 
the relative velocities of the reaction partners and of the 
target nuclear excited states states (par. (3) of Sect. [2Jl ■ In 
contrast, the provided rates are not corrected for stellar 
electron screening effects. 

3.1. Reaction rate calculations: general framework 

Some basics of the HF formalism adopted in the code 
MOST are briefly recalled in e.g. Arnould & Goriely 
(2003), and are not repeated here. The limits of its va- 
lidity are also reminded above. Let us just point out the 
following: 

(1) under local thermodynamic equilibrium conditions, 
the effective stellar rate of I + j — > L + k per pair of 
particles in the entrance channel at temperature T taking 
due account of the contributions of the various excited 
states /i of the target is expressed in a classical notation 
as 



ftWiC) = ( 



-) 

7rm/ 



1/2 



N A 



(fcT) 3 / 2 G(T) 



E 



E 



kT 



dE, (2) 



where k is the Boltzmann constant, m the reduced mass 
of the 1° + j system, Na the Avogadro number, a^ k {E) 
the cross section at relative energy E of the 1^ + j — > 
L + k reaction, and G(T) is the partition function given 
by Eq. □ where and are defined. The BRUSLIB 
photodisintegration rates are estimated by applying the 
reciprocity theorem on the radiative capture rates derived 
from Eq. [5] This procedure leads to 



(2Jj + l)(2J j + 1) Gi(T) /AjAj\ 3 / 2 



{2J L + 1) 

kT \ 3/2 



G L (T)\ Ai 



( kl V' AT / X 



U,i) 



-Qjy/kT 



(3) 



where Qj 7 is the Q-value of the I°(j,j)L° capture. Note 
that, in stellar conditions, the reaction rates for targets 
in thermal equilibrium obey reciprocity since the forward 
and reverse channels are symmetrical, in contrast to the 
situation which would be encountered for targets in their 
ground states only (Holmes et al. 1976). 

The uncertainties involved in any HF prediction are 
dominated by those involved in the evaluation of the nu- 
clear quantities necessary for the calculation of the cross 
sections, such as the masses, deformations, matter distri- 
butions, single-particle levels, and level densities of tar- 
get and residual nuclei, as well as the optical potentials. 
Special problems are also raised by the evaluation of the 
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photon widths T 7 . The interested reader is referred to e.g. 
Arnould & Goriely (2003) for more details. 

3.2. Nuclear masses, level densities, and partition 
functions 

The BRUSLIB MOST predictions rely on the experimen- 
tal nuclear mass data compiled by Audi et al. (2003). 
When not measured (a very common situation for vari- 
ous astrophysics applications), use is made of the HFB-9 
microscopic mass model (Goriely et al. 2005). This model 
also provides the necessary information on nuclear defor- 
mation, charge and matter distributions, pairing proper- 
ties and single-particle spectra. The nuclear level densi- 
ties are extracted from a microscopic model developed by 
Goriely (1996) (see also e.g. Demetriou & Goriely 2001). 

The nuclear partition functions Gj(T) entering the 
evaluation of the astrophysical reaction rates (Sects. [3]and 

are calculated from (i) a summation over experimen- 
tally known levels (Eq. ^) up to an excitation energy e w 
above which the knowledge of the energy spectrum is con- 
sidered to be incomplete, and (ii) a generalization of Eq.^ 
involving an integration over a level density evaluated as 
described above. 

3.3. Optical potentials 

Phcnomenological optical potentials (OPs) (generally of 
the Woods-Saxon type) may not be well suited for certain 
applications, particularly those involving exotic nuclei. 
It is considered profitable to use more microscopically- 
based potentials, whenever possible. A semi-microscopic 
OP, usually referred to as the JLM potential (Jeukenne et 
al. 1977), is available for the description of the nucleon- 
nucleus case. This OP has been revised recently for nu- 
cleons incident on spherical or quasi-spherical nuclei with 
masses 40 < A < 209 at energies ranging from 1 keV 
to 200 MeV (Bauge et al. 2001). The resulting new ver- 
sion gives a global satisfactory agreement with experi- 
mental data, even if some improvements would be most 
welcome, especially in the low-energy domain and in the 
treatment of deformed or exotic nuclei. It is adopted for 
the BRUSLIB rate evaluations. 

The situation for the a-particle-nucleus OPs is much 
less satisfactory, and one still has to rely on phenomeno- 
logical potentials. Most of the proposed OPs are derived 
from fits to elastic a-nucleus scattering data at energies 
E > 80 MeV or, in some cases, to (n,a) cross sections 
at lower energies. However, the OP, and in particular its 
imaginary component, is known to depend strongly on en- 
ergy below the Coulomb barrier. As a consequence, its ex- 
trapolation to sub-Coulomb energies that is necessitated 
by astrophysics applications is more insecure than in the 
case of nucleons. Several attempts to device a global a- 
nucleus OP for the description of the scattering and reac- 
tion cross sections at energies E < 20 MeV of better rele- 
vance to astrophysics have been conducted (e.g. Arnould 



& Goriely 2003 for more details and references). The 
scarcity of experimental data, particularly in the A > 100 
mass range, limits dramatically the predictive power of 
any of the constructed global OPs. This has the immediate 
consequence of reducing the reliability of the rate predic- 
tions as they depend sensitively on the a-particle-nucleus 
OPs. The BRUSLIB a-capture rates are calculated with 
the global OP III developed by Demetriou et al. (2002). 

3.4. j -ray strength function 

When applied to radiative captures, the total photon 
transmission coefficient entering the calculation of the cr^ 7 
cross section of Eq. [21 is dominated by the El transitions. 
The calculation of the iT-strength function necessitates 
the knowledge of the low-energy tail of the Giant Dipole 
Resonance (GDR) of the compound system formed in the 
reaction process. The photon transmission coefficient is 
most frequently described in the framework of the phe- 
nomcnological generalized Lorentzian model (e.g. Goriely 
1998). The Lorentzian GDR approach suffers, however, 
from shortcomings of various sorts. On the one hand, it 
is unable to predict the enhancement of the El strength 
at energies below the neutron separation energy demon- 
strated by nuclear resonance fluorescence experiments. 
This departure from a Lorentzian profile may manifest 
itself more clearly for neutron-rich nuclei, and especially 
in the form of a so-called pygmy El resonance (Govaert et 
al. 1998, Zilges et al. 2002). On the other hand, even if a 
Lorentzian function provides a suitable representation of 
the El strength, the location of its maximum and its width 
remain to be predicted from some underlying model for 
each nucleus. For astrophysical applications, these prop- 
erties have often been obtained from a droplet-type of 
model (Myers et al. 1977), which clearly lacks reliability 
when dealing with exotic nuclei. This introduces large un- 
certainties in certain rate estimates (e.g. Goriely & Khan 

2002, Goriely et al. 2004) 

In view of this situation, it is clearly of substantial 
interest to develop models of the microscopic type which 
are hoped to provide a reasonable reliability and predictive 
power for the ET-strength function. Various attempts of 
this sort have been conducted (e.g. Arnould & Goriely 

2003, and references therein). The BRUSLIB MOST rate 
calculations are based on the semi-microscopic QRPA El 
model developed by Goriely et al. (2004). 

4. A confrontation between measured and 
calculated reaction rates 

In order to evaluate the overall quality of the BRUSLIB 
reaction rate predictions, this section presents a confronta- 
tion between selected experimental data and calculations 
in which only the ground state contribution (fi = 0) is 
taken into account in Eqs. and |3J the consideration of 
target excited states being irrelevant in the laboratory 
conditions. 
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Fig. 1. Comparison of BRUSLIB Maxwellian- averaged 
(11,7) rates (eru)th with experimental values (Bao et al. 
2000) at T = 3.5 x 10 8 K 



Figure ^ compares the BRUSLIB predictions of the 
Maxwellian-averaged (n,7) rates (av) at T — 3.5 x 10 8 K 
with experimental data for some 228 nuclei heavier than 
40 Ca included in the compilation of Bao et al. (2000). It 
appears that the calculations agree with all data to within 
a factor of three. Figure[3]compares the experimental cross 
sections for some low-energy (p,7) reactions on targets 
heavier than Fe with the corresponding BRUSLIB rates. 
Again, and broadly speaking, the agreement is very sat- 
isfactory for the vast majority of the (p,7) data. It is dif- 
ficult, however, to be much more specific, as the quality 
of this agreement may depend more or less drastically on 
temperature, in contrast to the situation encountered in 
the neutron capture case. In addition, the comparison is 
limited to targets up to Sn. Experiments with heavier nu- 
clei would be most welcome, as they might help refining 
the predictions. 

The situation is by far less clear for the (a, 7) reac- 
tions. This results from the lack of a large enough body 
of experimental data for sub-Coulomb cross sections com- 
bined with the difficulties to construct global and reliable 
cv-nucleus OPs fSect. l3~3^) . These theoretical problems are 
magnified by the fact that, at the sub-Coulomb energies 
of astrophysical relevance, the reaction rate predictions 
are highly sensitive to these potentials through the cor- 
responding cv-particle transmission coefficients. Figure [3] 
displays a comparison between some low-energy measure- 
ments of (a, 7) cross sections and MOST predictions used 
for evaluating the BRUSLIB rates. The quality of the 
agreement appears to vary from case to case, and is also 
highly sensitive to temperature. Uncertainties in the cal- 
culated values originate mainly from those in the nuclear 
level densities and a-nucleus OPs (see Arnould & Goriely 
2003 for details, and in particular for a discussion of the 
144 Sm (a , 7) 148 Gd reaction of special astrophysics inter- 
est). 




1 1.5 2 2.5 3 3.5 4 1 1.5 2 2.5 3 3.5 4 



E[MeV] E[MeV] 



Fig. 2. Comparison between the BRUSLIB predictions 
and the measured cross sections of some (p,7) reactions 
on targets heavier than 84 Sr (see Arnould & Goriely 2003 
for references to the laboratory work) 

5. The nuclear network generator NETGEN 

The Nuclear Network Generator NETGEN is an interac- 
tive Web-based tool to generate Maxwellian-averaged nu- 
clear reaction rates for networks and temperature grids 
specified by the user. It is fully documented at the web 
address http://www-astro.ulb.ac.be/Netgen 

NETGEN relies mainly on the BRUSLIB /NACRE li- 
brary. It also makes use of 

(1) the post-NACRE compilations by Iliadis et al. (2001) 
and Descouvemont et al. (2004) (see Table 1); 

(2) experimentally-based published rates for about 
70 charged-particle induced reactions not included in 
NACRE. Specific references are provided for each of these 
reactions; 

(2) more than 200 experimentally-based radiative neutron 
capture rates. Most of these are adopted from Bao et al. 
(2000). Specific references are given in each case; 

(3) /3-decay and electron-capture rates, including (i) labo- 
ratory measurements compiled by Horiguchi et al. (1996), 
(ii) theoretical estimates based on the evaluation of indi- 
vidual transitions (Takahashi & Yokoi 1987), on the gross 
theory (Tachibana et al. 1990), or on the more microscopic 
ETFSI + cQRPA model (Borzov & Goriely 2000). 
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E cJ MeV > E cm( M eV) E cm (MeV) 

Fig. 3. Cross sections for the (a, 7) reactions on 37 C1, 
42 Ca, 56 Fe, 58 Ni, 62 - 64 Ni, 70 Ge, 106 Cd and 144 Sm (see 
Arnould & Goriely 2003 for references to the laboratory 
work) . Black and open circles represent experimental data. 
The solid lines give the calculated cross sections used to 
construct the BRUSLIB rates 

For the sake of completeness, the references of Table 1 
that are not contained in the NACRE, Iliadis et al. (2001) 
or Descouvemont et al. (2004) compilations are also in- 
cluded in NETGEN. 

For each reaction or /3-decay rate, NETGEN selects 
by default the data source that is considered to be the 
most reliable. We select in order of preference the lat- 
est available compilation (NACRE, Iliadis et al. 2001, 
or Descouvemont et al. 2004), experimental data, de- 
tailed microscopic calculations, and BRUSLIB rates de- 
rived from global calculations. The user may nevertheless 
adopt another choice for selected cases by specifying 'bibli- 
ographic indexes' for each reaction, the table of rates being 
accompanied by a 'log file' listing the selected data source 
among all those available in the library. A FORTRAN 
program handling the rates is also made available. Note 
that all rates duly take into account the contribution from 
the excited states of the target nuclei, as discussed in the 
previous sections. 

Various NETGEN options are currently offered 
through the web interface, as described in detail at 
http: //www- astro, ulb.ac. be/Netgen 

(1) Generate a table of reaction rates on a temperature 
grid for a network that has been 

— typed in reaction by reaction (offering the possibility 
to select non-default rates, as mentioned above), or 

— generated automatically between interactively-selected 
boundaries on the proton and mass numbers, and in- 
volving various possible sets of reactions (p-, n-, a- 
captures, /3-decays and/or photodisintegrations), or 



— uploaded on the server. 

(2) Plot individual reaction rates, and provide .gif or .ps 
files. 

6. Conclusions 

This is the first release in an astronomy and astrophysics 
journal of the BRUSLIB nuclear reaction rate library and 
of the nuclear network generator NETGEN. The format 
of the packages is chosen in order to make them easily 
accessible to astrophysicists and to be well suited for a 
large variety of astrophysics needs. They are made avail- 
able through the web site http://www-astro.ulb.ac.be 

The BRUSLIB NACRE package contains a detailed 
experimentally-based evaluation and compilation of the 
rates of 86 proton or a-capture reactions on (mainly) sta- 
ble targets up to Si for temperatures ranging from 10 6 to 
10 10 K. The electronic files contain much more informa- 
tion than published by Angulo et al. (1999). 

The NACRE data are complemented with about 
100000 thermonuclear rates of nucleon and a-captures on 
about 8000 8 < Z < 110 nuclei located between the pro- 
ton and neutron drip lines. The calculations are based 
on a statistical Hauser-Feshbach model featuring a micro- 
scopic (or at least very close to microscopic) evaluation of 
the basic ingredients of the model. These predictions are 
seen to compare favourably with the limited set of exper- 
imental reaction cross section data on intermediate-mass 
and heavy nuclei at energies close to those of astrophysics 
relevance. The rates of photodisintegration of the whole 
set of nuclei are also provided. They are derived from the 
application of the reciprocity theorem. 

NETGEN is an interactive web-based tool allow- 
ing the construction on a user-friendly basis of nu- 
clear reaction networks specified by the user on a 
temperature grid of his/her choice. A full documen- 
tation of its use can be found at the web address 
http://www-astro.ulb. ac. be/Netgen 

It is hoped that the easy availability of a very large 
set of nuclear reaction rate evaluations and predictions, 
as well as of a nuclear network generator will be helful to 
many researchers for a large variety of applications. Quite 
clearly, this enterprise is of a highly dynamical and long- 
term nature. BRUSLIB will be continuously improved and 
expanded, and new releases will be made accessible to the 
community every time a substantial enough body of new 
data becomes available. 
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